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Abstract. Using detailed exact results on pair-correlation functions of Z-invariant 
Ising models, we can write and run algorithms of polynomial complexity to obtain 
wavevector-dependent susceptibilities for a variety of Ising systems. Reviewing recent 
^ ' work we compare various periodic and quasiperiodic models, where the couplings 

and/or the lattice may be aperiodic, and where the Ising couplings may be either 
I ferromagnetic, or antifcrromagnetic, or of mixed sign. We present some of our results 

' on the square-lattice fully-frustrated Ising model. Finally, we make a few remarks 

on our recent works on the pentagrid Ising model and on overlapping unit cells in 
three dimensions and how these works can be utilized once more detailed results for 



o 

\o 
o 

' pair correlations in, e.g., the eight-vertex model or the chiral Potts model or even 

three-dimensional Yang-Baxter integrable models become available. 

■ PACS numbers: 05.50.+q, 64.60.Cn, 75.10.-b 



■ 1. Introduction 

Detailed exact results on pair-correlation functions and susceptibilities are usually hard 
to come by, even in models that are "exactly solvable." The best results are available for 
Z-invariant Ising-type models, and even though the results are often not in closed form, 
we may follow Tony Guttmann and claim an exact result if we can produce an algorithm 
of polynomial complexity to obtain detailed information on, say, wavevector-dependent 
susceptibilities for a variety of Ising systems. 

For a long time, the paper of Wu, McCoy, Tracy and Barouch P] has been the 
standard work on the pair correlations and susceptibility of the square-lattice Ising 
model. Only in the last few years do we know how to do much better. We now have 
algorithms to construct long high- and low-temperature series together with detailed 
expansions at criticality for the susceptibility O Ej using quadratic recurrence relations 
[4j for pair correlations. We also have much more analytic and numerical information for 
pair correlations and susceptibilities in more general Z-invariant inhomogeneous Ising 
models [SIISIIZ!, together with new understanding from a field theory approach ^j. 
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1.1. Baxter's Z-invariant inhomogeneous Ising model 

We can start with an inhomogeneous Ising model on a square lattice with reduced 
interaction energy 

(31~L ^ ^ (-^m,n'^m,n'^m, 71+1 ~l~ ■^m,n^m,n^m+l,n) (l) 



where (3 = l/k^T the inverse temperature. At site {m,n) the spin takes values 
am,n = ±1, while K = (3J and K = (3J are "horizontal" and "vertical" dimensionless 
coupling constants. 

One important quantity to study is the wavevector-dependent susceptibility x(ct) 
defined by 

ii^Z^E^'^'^^'^^K^r^rO - (2) 



fcBTx(q) = x(q) 



This is the Fourier transform of the connected pair correlation function. In >C is the 
number of lattice sites, r and r' run through all sites and q = {qx,Qy) is the wavevector. 

Baxter's Z-invariant Ising model jHl CSl is defined in terms of oriented rapidity 
lines on which the rapidity variables live. (The rapidity lines can be moved at will, 
without changing the partition function Z, so we do not have to take a square lattice.) 
The faces bounded by the rapidity lines are alternatingly colored black or white and 
spins live somewhere on the black faces forming the Ising system {cTj.}, whereas dual 
spins live on the white faces forming the dual Ising system {a**}. Spins interact only 
with nearest neighbours from which they are separated by precisely two rapidity lines, 
see Fig. ^ The couplings K and K are parameterized in terms of elliptic functions of 
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Figure 1. Two oriented rapidity lines with rapidities Ui and Vj are drawn dashed. 
The orientation of the couphng with respect to these rapidity hnes defines whether we 
have (a) horizontal coupling Kij ~ K{ui,Vj) or (b) vertical coupling Kij — K{ui,Vj). 



modulus k, i.e. 



sinh {2K{ui, U2)) = k sc(-ui — U2, k') = cs(K(A;') + U2 — ui, k') 



sinh {2K{ui, U2)) = cs(mi — U2, k') = k sc(K(A;') + U2 — Ui, k') (3) 
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where 

k' = Vl — k'^, sc{v,k) = sn{v,k)/cn{v,k) = l/cs{v,k) (4) 

and K(/c) is the complete elhptic integral of the first kind. 

We note that K and K are interchanged if we replace Ui by U2 ± K{k') and U2 by 
Ui. In other words, flipping the orientation of a rapidity line j is equivalent to changing 
its rapidity variable uj to uj ± K{k'). 

1.2. Two-point correlation functions 

Z-invariance implies jH] that the pair correlation only depends on elliptic modulus k 
and the values of the 2m rapidity variables mi, . . . , U2m that pass between the two spins, 
implying the existence of an infinite set of universal functions g2,g4:, ■ ■ ■ ,g2m, ■ ■ ■ such 
that for any permutation P and rapidity shift v 

{aa') = g2m{k; Mi, • • • , U2m) = g2m{k; Up(i) + . . . , Mp(2m) + v). (5) 

Here Uj = Uj if the jth rapidity line passes between the two spins a and a' in a given 
direction and uj = Uj + K{k') if it passes in the opposite direction [TT]. 

If two of the rapidity variables passing between the two spins differ by K{k'), they 
can be viewed as belonging to a single rapidity line moving back and forth between 
these two spins, i.e. 0, 

g2m+2{k;Ui, . . . ,U2m,U2m+l,U2m+l + K(A;')) = g2m{k;Ui, . . . , U2m) ■ (6) 

1.3. Quadratic identity for pair correlation 

For general planar Ising models we can derive a quadratic relation between pair 
correlation functions 4j. For the situation in Fig. |21we have 




^1 

Figure 2. At crossing 1 of two rapidity lines we have spin sites xi and yi together 
with dual spin sites x| and y^. At crossing 2 we have positions X2 and y2 together 
with and y^ oriented similarly. 

sinh(2iri) sinh(2ir2){(o"xiO"x2)(o"yiO"y2) - {a^^cTy.^) {(yy.a^^)} 

+ {(^xi^x^)*KifXy*)* - {cr^iayr,yWy,a,,y} = 0, (7) 
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with two arbitrary nearest-neighbour pairs of spins at the sites {xi,yi} 7^ {x2,j2}, and 
corresponding nearest-neighbour pairs of dual spins at {x^,y^} and {x2,y2}, whereas 
sinh(2i^^j) sinh(2ii'*) = 1, {i = 1, 2). Orientations have to be consistent as in the picture, 
otherwise the plus changes to a minus. 

Restricted to Z-invariant Ising models, the quadratic identity reduces to 

A;^sc(m2 — Ml, k')sc{u4 — U3, k') 

X {g{ui, U2, Us, U4, ■ ■ ■) g{- ■ ■) - g{ui, U2, ■ ■ •) g{u3, M4, ■ ■ ■)} 

+ {g*iui,U3, ■ ■ ■) g*{u2, M4, ■ ■ ■) - 9*iui, Ui, ■ ■ •) g*{u2, M3, ■ ■ ■)} = 0, (8) 

with "■ ■ ■" short-hand for all other rapidity variables Ms, Me, ■ ■ common to all g^s and 
g*'s (passing between all eight sites). 

Knowing g{u, u, ■ ■ ■ ,u) and g{v, u, ■ ■ ■ ,u), with all or all-but-one of the rapidities 
equal, all other g^s can be calculated by recurrence. Therefore, the knowledge of the 
diagonal and next-to-diagonal pair correlations in the uniform asymmetric {K 7^ K) 
square- lattice Ising model suffices jHl Ej . 

Jin has found the scaling form for the general Z-invariant case in the critical regime 
fsi 1 [3 EI ■ From this we can find the first two terms of the susceptibility in several 
lattices confirming Guttmann's extended lattice-lattice scaling jT2l HH] ■ 

1.4- Summary of findings 

We can make the couplings J and/or the lattice aperiodic. We find the following 
001111113]: 

• Periodic lattice with periodic couplings: Periodic xi^)^ with peaks at sites 
commensurate with reciprocal lattice, becoming sharper and sharper as T ^ Tc. 
This includes fully-frustrated cases. 

• Periodic lattice with ferromagnetic couplings varying quasiperiodically: Periodic 
x(q), with peaks at reciprocal lattice sites, sharper and sharper as T — > Tc 

• Periodic lattice with mixed FM and AFM couplings quasiperiodically arranged: 
Periodic x(ct), with more and more incommensurate peaks within unit cell as 

• Quasiperiodic lattice: Quasiperiodic x(ct), with more and more peaks visible closer 
to Tc 

For Z-invariant lattices, we can evaluate numerically to high accuracy from 
the recurrence relations for the pair correlations. However, the structure is clearer in 
density plots E] . 

1.5. Generalized Fibonacci Ising lattices 

We can assign the couplings according to quasiperiodic sequences in horizontal, 
vertical, and/or diagonal directions. We can use ^1] de Bruijn's generalized Fibonacci 



Pair Correlations and Susceptibilities in Exactly Solvable Models 



5 



sequences, assigning different couplings according to the sequence of zeros and ones 

Pj{n) = [7 + + l)/aiJ - L7 + ^/«iJ (9) 

with 

«. = ^[(J + 1) + V(J + 1)^ + 4] (10) 

We found that such ferromagnetic cases differ very httle from periodic cases. 

For the mixed ferro/antiferro case, the depends strongly on the sequence 

chosen, even for the simplest examples just adding signs to the couplings of the Onsager 
square lattice model by gauge transform . 

1.6. Pentagrid Ising lattices 

Following Korepin, we can use de Bruijn's pentagrid for the rapidity lines leading to a 
spin model and its dual taking alternating sites of a Penrose tihng [El- By the quadratic 
recurrence relations we can compute a big collection of pair correlations. Then, in order 
to calculate we have developed a novel way of determining the pair probability 

of local environments on a Penrose tiling, which can also be used once more detailed 
results for pair correlations in e.g. the eight-vertex model or the chiral Potts model 
become available. Full details will be published elsewhere |15j . 

2. Fully- frustrated square-lattice Ising model 

There are two common ways to define a fully-frustrated Ising model on a square lattice. 
One way is the checkerboard version of coloring the faces of the lattice alternatingly 
black and white and taking three of the couplings around each black square +J and the 
fourth one —J. The other way is to take all horizontal couplings +J and the couplings 
in vertical columns alternatingly + J and — J. This is illustrated in Fig. El 



+J -J 



Figure 3. Two versions of the fully-frustrated square-lattice Ising model with 
ferromagnetic couplings +J and antiferromagnetic couplings — J: (a) Checkerboard 
version on the left; (b) version periodic in vertical direction on the right. 
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It must be noted that these two versions do not define fundamentally different 
models. They are, in fact, related with each other by a simple gauge transformation 
of fiipping the signs of each second horizontal pair of rows of spins, as is indicated in 
Fig. m Of course, there are many other ways of fiipping subsets of the spins and one 



Figure 4. Gauge equivalence of the checkerboard and vertically periodic cases: Flip 
the signs of the red spins to go from the one to the other, or back. 

can define infinitely many other related less regular frustrated models. 

There are several approaches to the explicit evaluation of pair correlations and 
susceptibilities in this model. 

2.1. Sum out every other spin 

One approach is to first sum out every other spin in order to arrive at an effective Baxter 
eight- vertex model [IHIIIZI, with diagonal couplings J and J' and four-spin coupling J4. 
If one does this for the checkerboard case, see also Fig. El one finds 



o 



-t>. 



a 



Figure 5. Spin cto to be "decimated out" with its four neighbours. 



1 



3^^' = + 1, e' 



g2 



(11) 



V2S^ + l' ' V2S^ + 1 

with K = J /k-QT, K' = J'/k^T and K4 = J^i/k-QT, while using the short-hand notation 



S = sinh 



2J 



sinh 2K 



(12) 
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In terms of Baxter's a, 6, c, d, this means 



a = b= ^/S^Tl, c=l, d = V2S^ + 1 (13) 

Note that only diagonal and four- spin interactions are left. All other interactions cancel 
out. 

This determines the mapping of Boltzmann weights. For correlation functions we 
have also to find the mapping of the spins involved. Spins ai, o^-, C3, in Fig.Elare left 
unchanged. However, the central spin oq in Fig. El has to be replaced after "decimation 
summation" according to 

5' ((71 + 0-3 + 0-2-0-4) A S'^(l - 0-iO-30-20-4)\ 1^.. 

"° 2^S^ + 1 2(2^^ + 1) J 

2.2. Partial duality transform 

After this, we first flip the sign of each 0-2O-4, interchanging a d, b c. We then have 

e'^ = e'^' = V2S^, e'^^ = ^^^^, ^ - sinh ^ (15) 

Next, we apply a partial Kramers- Wannier duality transform on every other spin 
[T6| IT7j . This leads to an Ashkin-Teller model, with new spins r on the same positions 
as a quarter of the original spins a and three quarters of the original positions empty. 
All original spins a are now expressed in terms of a's and r*'s (disorder variables, equal 
to other original a's). 

Finally, replacing a — > ar, r — > r, a* — > a*, r* — > a*T*, the Ashkin-Teller 
model factorizes into two Ising models, at dual temperatures = Pa Jo, Kr = 13* Jq, 
with 

sinh(2ir<,) = = . , , = Vk (16) 

^ ^ 52 + 1 + V252 + 1 smh(2ir,) ^ ^ 

2.3. Difference equations for the correlation functions 

Writing C{m,n) = {crQfi(^m,n) C{m,n) = {TQflTm,n) Xr -^^^ resulting 
independent Ising models, we can determine these recursively |31 [7j from 

[C{m,n + l)C{m,n-l) - C{m,nf] 

+ k[C{m + l,n)C{m-l,n) - C{m,nf] = (17) 
[C{m + l,n)C{m-l,n) - C(m,?2)^] 

+ k[C{m,n + l)C{m,n-l) -C{m,nY] = (18) 
[C{m,n)C{m+l,n+l) - C(m+ 1, n)C(m, n+1)] = 

k [C{m,n)C{m + l,n+l) - C{m+l,n)C{m,n + l)] (19) 
Vk [C(m+1, n)C{m — l, n) + C(m — 1, n)C'(m + l, n) 

+ C{m,n + l)C{m,n-l) + C{m,n-l)C{m,n + l)] 

= {k + l)C{m,n)C{m,n) (20) 
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C(0,0) =C(0,0) = 1, C(l,0) + v^C(0,l) = v/A^ (21) 
C(m,n) = C{n,m) = C{\m\, \n\) 

C{m, n) = C{n, m) = C{\m\, \n\) (22) 
2.4- Pair correlation function of fully-frustrated model 

Finally, the pair correlation functions of the original model of Fig. |21(a) are given by 

i^ki '^k+2m,i+2n) = ("1)" Cim, u) C(m, u) (23) 

(^fe«^fe+2m-l,«+2n-l) = (24) 
i^kl ^ k+2m-l,l+2nl 

-.{C{m-l,n)C{m,n) + C{m,n)C{m-l,n)) (25) 



2V2S2 + 1 

i^kl ^k+2m,l+2n~l) 
±{-lYS 



(C(m, n-l) C{m, n) + C(m, n) C{m, n - 1)) (26) 

In the last equation we have to choose plus (+) if min(A; + /, /c + 2m + I + 2n — \) = 
even and minus (— ) if mm{k + l,k + 2m + / + 2n — 1) = odd. Therefore, its contribution 
to x{Qx, %) averages out to zero. If we had started from the vertically periodic case 
Fig. 121(b), we must omit all three (—1)"" factors in the above. 

Obviously both models have a periodic x{Qxi %) with commensurate peaks only, as 
can be illustrated with density plots. Using the results of this subsection it is also not 
difficult to derive a long series expansion along the lines of the Ising model work 0. 



2.5. Other approaches 

The square-lattice fully- frustrated Ising model is dual to the Ising model in field mk-oT /2 
introduced by Lee and Yang. For this case it in convenient to view a square of four 
spins as a vertex of a 16- vertex model and then to apply gauge transforms llj. Now 
the model decouples as a product of two Ising models at the same temperature. For 
the special case of the square-lattice dimer model these two models are at the critical 
temperature and the explicit formula for the two-point function has been published ^H] 
and the monomer-monomer correlation function and its Fourier transform x{Qxi %) have 
been studied in some detail by Kong pTJj . 



2. 6. Final remarks: Checkerboard Ising and chiral Potts 

The physical free-fermion model corresponds to a checkerboard Ising model ^U] with 
real or imaginary elliptic modulus k. In general, the correlation functions satisfy 
a 2-by-2 matrix generalization [201 difference equations used above, allowing 

incommensurate solutions. 
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For the two-dimensional integrable iV-state chiral Potts model a similar mapping 
to an A^-state generalization of the free-fermion 8- vertex model exists This model 
is a most natural generalization of Onsager's two-dimensional Ising model to more than 
two states per spin. The pair-correlation functions have not yet been solved for this 
model; only the conjecture for its order parameters (one-point functions) has been proved 
recently by Baxter [221 ■ The model is on a submanifold in the commensurate phase of the 
more general chiral Potts model. The full chiral Potts model also has incommensurate 
phases. However, Jin's results [SB] do not seem to support the presence of a Lifshitz 
point in the classical two-dimensional model. 

3. Overlapping unit cells 

Gummelt [24 motivated by physical considerations, has proposed a description of 
quasicrystals in terms of overlapping unit cells, with the regular Penrose tiling described 
by overlappings of decorated decagons. We have used a multigrid method based on 
de Bruijn's work to produce a new example of 3-dimensional overlapping unit cells, 
quasiperiodic in two directions and periodic in the third. Full details are presented 
elsewhere |25j . 

As our construction is based on a multigrid of five grids of parallel planes, resulting 
from a projection of the five- dimensional hypercubic lattice into three dimensions, one 
may construct aperiodic integrable three-dimensional models with spectral variables 
living on the grid planes and with Boltzmann weights satisfying the tetrahedron 
equations, generalizing the pentagrid Ising model construction. 
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